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SIGNIFICANCE  AND  EXPLANATION 

Many  boundary  value  problems  for  elliptic  equations  can  be  formulated  as 
variational  problems  in  which  a quadratic  functional  must  be  minimized  on  the 
subspace  of  functions  which  satisfy  the  boundary  conditions.  The  best  known 
example  of  such  variational  problems  arises  when  Laplace's  equation  must  be 
solved  subject  to  prescribed  boundary  values,  in  which  case  the  variational 
problem  requires  the  minimization  of  the  Dirichlet  integral.  This  provides 
the  starting  point  for  the  widely-used  finite  element  method  for  computing 
solutions  of  potential  problems  such  as  occur  in  steady  state  heat  conduction 
or  inviscid  fluid  flow,  for  example. 

An  elliptic  variational  inequality  is  a generalization  of  variational 
problems  of  the  above  type:  it  is  required  to  minimize  a quadratic  functional 
on  a convex  set.  Elliptic  variational  inequalities  arise  in  contexts  in  which 
a physical  system  is  subject  to  restraints,  as  when  a membrane  is  stretched 
over  an  obstacle.  (If  no  obstacle  is  present,  the  deflection  of  the  membrane 
is  governed  by  Laplace's  equation.  The  obstacle  is  the  new  feature  that  leads 
to  the  variational  inequality.  The  region  of  contact  between  the  membrane  and 
the  obstacle  is  not  known  beforehand.) 

This  report  first  describes  the  theory  of  approximating  elliptic  varia- 
tional inequalities,  and  then  goes  on  to  discuss  the  numerical  solution  of 
such  problems  using  finite  elements. 
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FINITE  ELEMENTS  AND 


VARIATIONAL  INEQUALITIES 


R.  GLOWINSKI* 


INTRODUCTION 

During  these  last  years  many  works  concerned  with  the  Numerical  Analysis  of 
Variational  Inequality  problems  have  been  published.  We  would  like  in  this 
paper  to  describe  some  of  the  results  which  have  been  obtained  with  emphasis 
on  the  Finite  Element  Approximation  of  these  Inequalities. 

The  content  is  as  follows  : 

1.  Abstract  Elliptic  Variational  Inequalities.  Existence,  Uniqueness, 
Approximation. 

2.  Specific  examples  and  error  estimates  for  conforming  finite  element 
methods . 

3.  Iterative  methods  for  solving  variational  inequalities. 

4.  Approximation  of  the  obstacle  problem  by  mixed  finite  element  methods. 

5.  Further  comments.  Conclusion. 
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I . ABSTRACT  ELLIPTIC  VARIATIONAL  INEQUALITIES.  EXISTENCE,  UNIQUENESS , 
APPROXIMATION. 

1.1.  Orientation. 

In  Sec.  1 we  just  consider,  following  GLOWINSKI  [53,  Ch.  1,2,3]  some  simple 
classes  of  Elliptic  Variational  Inequalities  (E.V.I.)  and  their  approximation 
by  Galerkin  type  methods.  More  specific  examples  will  be  considered  in  Sec.  2. 

1.2.  Functional  context. 

We  introduce  the  following  mathematical  objects  : 

. A real  Hilbert  space  V equipped  with  the  inner  product  (•,♦)  and  the  corres- 
ponding norm  1 1 • ||  . 

. A bilinear  continuous  form  a : V*V  -+  R,  V-elliptic  (i.e.  3a  > 0 such  that 
2 

a(v,v)  >a||v||  Vv  e V)  ; we  don’t  assume  that  a(*,*)  is  symmetric. 

. A form  L : V ■*  R,  linear  and  continuous. 

. A closed,  convex,  non  empty  subset  K of  V. 

. A functional  j : V -*■  R (=  R u{+«>}  u { — <»} ) , convex,  lower  semi  continuous 
(l.s.c.),  proper  (i.e.  j(v)>-°°  VveV,  j £ +<*>) . 

1.3.  Two  classes  of  E.V.I. 

Let  us  consider  now 

Find  ue  K such  that 

(EVI)j 

a(u,v-u)  > L(v-u)  V v e K, 

and 

Find  ueV  such  that 

(evi)2 

a(u , v-u)  + j (v) -j  (u)  > L(v~u)  VveV. 

In  the  sequel  (EVI) j (resp  (EVI)2)  will  be  denoted  as  an  EVI  of  the  first 
(resp.  second)  kind. 

Remark  1.1.  : We  can  find  in  the  litterature  more  complicated  EVI,  (cf.  fv 
example  LIONS  [65])  and  also  the  generalization  called  Quasi  Variational  In 
equalities  (QVI)  introduced  by  BF.NSOUSSAN-LIONS  (see  for  example  BF.NSOUSSAN- 
L10NS  [ 9 ] , BAIOCCHI-CAPELO  T 6]  and  the  bibliography  therein). 
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Remark  I . 2 : If  K=V  in  (EVI) ^ and  j - 0 in  (EVI)^,  then  both  problems  reduce 
to  the  standard  linear  variational  equation 


(1.1) 


Find  u e V such  that 
a(j,v)  = L(v) V v e V. 


Remark  1 . 3 : The  distinction  between  (EVI) ^ and  (EVI^  is  rather  artificial 
(theoretically  at  least)  since  (EVI)  is  equivalent  to 

(Find  ue  V such  that 

a(u,v-u)  + IK(v)-IR(u)  ^ L(v-u)  Vv  e V 

where  1^  is  defined  by 

!0  _if  v e K, 

+“  if  v { K. 


The  functional  IR  is  called  the  indicator  functional  of  K,  and  since  K is  a 
closed , convex , non  emp ty  subset  of  V,  1^  is  l.s.c.,  convex,  proper. 
Therefore  (EVI)  { is  a special  case  of  (EVI),,  ; however  formulation  (EVI) 
is  usually  more  practical. 


Remark  1 . 4 : If  a(*  ,*)  is  symmetric  (EVI) | and  (EVI^  are  respectively 
equivalent  to  the  minimization  problems 


( "IT  | ) 
where 


Find  u e K such  that 
J(u)  ^ J(v)  V v e K 


(1.3) 


J(v) 


1 

2 


a(v,u)-L(v) , 


and 


<v 


Find  u e K such  that 
J(u)+j(u)  ^ J(v)+j(v) 
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where  J(*)  is  still  defined  by  (1.3). 


1.4.  Existence  and  Uniqueness  results  for  (EVI) j , (EVI^. 

From  LIONS-STAMPACCHIA  [66]  we  have  the  following 

Theorem  1.1  : If  the  above  hypotheses  on  V,a,L,K  hold  then  (EVI) j has  a 
unique  solution. 

Proof  : (1)  Uniqueness 


Let  U| 

and  u2  be  two  solutions. 

Then 

(1.4) 

a(Uj ,v-u()  > L(v-Uj) 

VveK,  Uj  e K, 

(1.5) 

a(u2,v-u2)  > L(v-u2) 

Vv  e K,  u2  e K. 

Taking  v=u2  in  (1.4),  v=Uj  in  (1.5)  we  obtain  by  addition  and  using  the 
V-ellipticity  of  a(*,*)  that 

a II  u2-u]  II  2 * a(u2-Uj  .u2-u,)  ^ 0 

which  implies  the  uniqueness. 

(2)  Existence 

It  is  known  from  the  Riesz  representation  theorem  that  there  exists  Ae  £ (V,V) 
and  Z e V such  that 

(1.6)  a(u,v)  = (Au,v)  Vu.veV, 

(1.7)  L(v)  = (£,v)  Vv^V. 


Then  if  u is  solution  of  (EVI) j we  have 
( (Au,v-u)  ^ (£,v-u)  VveK, 


(1.8) 


| u e K 

which  is  equivalent  to  (EVI) ^ . Then  (1.8)  is  equivalent  to 


4 

Li 
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(1  .9) 


(u-p  ( Au-2.)  -u  ,v-u)  <0  Vve  K, 


| u e K , p > 0 , 
and  (1.9)  is  equivalent  to 

(1.10)  u = PK(u-p(Au-Jl)  , p>0  , 

where  in  (1.10)  PK  is  the  projection  operator  from  V to  K in  the  ||  • || -norm. 

It  follows  from  (1.10)  that  every  solution  of  (EVI) j is  also  solution  of 
the  fixed-point  problem  (1.10)  for  any  p > 0,  and  conversely  if  there  exists 
a particular  p such  that  (1.10)  has  a solution,  this  solution  is  also  solu- 
tion of  (EVI) j . 

A sufficient  condition  for  (1.10)  to  have  a solution  is  that  the  mapping 
from  V to  V 

v -*•  Pk(v-p(av-H)) 

is  strictly  and  uniformly  contracting  for  p well  chosen. 

Let  denote  by  w the  vector  P„(v-p (Av-£) ) and  let  consider 

K. 

Wi  = PR(vi-p(Avi-£))  , i=l ,2  . 

Since  the  projection  mapping  P^  is  a contraction,  we  have 

(1.11)  l|w2-wj|  5 ||v2-v)-pA(v2-v1)||  . 

From  (1.6),  ( 1 . I 1 ) we  obtain 

||  w2_w|  II  2 5 If  v2_v|  H 2 “2  a(v2"vi  »v2~vi)+p2  ||  A.Cv2-v,  ) II  2 ^ 

- li  v2_v,  ||  2 -2pct  ||  v2-Vl  ||  2 -*-p2  H A||  2 ||  v2-v,  ||  2 

so  that 

(1.12)  liw2-w,  || 2 S (l-2pa  + p2  ||A||2  ) ||v2-Vl||2  . 
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If  p obeys  (1.13),  the  fixed  point  problem  (1.10)  has  a solution  which  is 
also  the  solution  of  (EVI),  and  we  know  that  this  solution  is  unique . 

Remark  1 . 5 : The  proof  of  Theorem  1.1  suggests  the  following  algorithm  for 
solving  (EVI) j : 

(1.14)  u°  e V,  given 
then  for  n > 0 

(1.15)  un+1  = PK(un-p(Aun-H) ) . 

2ot 

From  the  proof  of  Theorem  1.1  it  follows  that  for  0<p  < r,  the  sequence 

II  » II  2 

n . 11 " 

(u  )n  defined  by  ( 1 . 1 4) , ( 1 . 1 5)  converges  strongly  in  V to  the  solution  u of 

(EVI),  . 

Practically  the  interest  in  ( 1 . 1 4) , ( 1 . 1 5)  is  quite  limited  in  most  applications 
(at  least  in  the  above  form)  since  we  usually  don't  know  Z or  A,  and  to  project 
on  K is  in  most  cases  a very  complicated  operation.  We  should  observe  that  if 
a(  • , • ) is  symmetric , then  A is  also  symmetric  and  (1.14), (1.15)  is  a gradient 
with  projection  algorithm  ; cf.,  e.g.,  CEA  [28!  for  a study  of  these  methods.  m 

Concerning  (EVI^  it  follows  from  LIONS-STAMPACCHIA,  loc.cit.,  that 

Theorem  1 . 2 : If  the  above  hypotheses  on  V,a,L,j  hold  then  (EVI)?  has  a unique 
solution. 

We  refer  for  the  proof  to  LIONS-STAMPACCHIA,  loc.  cit.,  and  also  to  GLOWINSki 
[53]  ; in  fact  in  this  proof  which  is  a variant  of  the  proof  of  Theorem  1.1, 
one  still  uses  a fixed  point  technique. 
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1.5.  Internal  Approximation  of  ( EV I ) j 

The  assumptions  on  V,a,K,L  are  those  of  Sec.  1.2. 


I 


- 


1.5.1.  Approximation  of  V and  K. 

The  parameter  h converging  to  zero,  the  space  V is  approximated  by  the  family 
(Vh)h  where  the  V^  are  closed  subspaces  of  V (usually  the  are  finite- 
dimensional) . 

Then  the  convex  set  K is  approximated  by  (K^)^  where  the  are  closed  convex 

subsets  of  V.  ; we  do  not  assume  that  K,  cR,  We  do  assume  however  that 
h h 

has  the  two  following  properties 


(i)  If_  (v^)^  is  such  that  v^  e Vh  then  all  the  weak  cluster  points  of 
(vh)h  belong  to  K, 

(ii)  There  exist  x/x  = K and  rh  : ^ *Si  such  that 

lim  rvv  = v strongly  in  V.Vvey 
h-*-0  “ 

Remark  1 . 6 : If  K^cK  Vh  , then  (i)  is  automatically  satisfied. 


1.5.2.  Approximation  of  (EVI) j . 
We  approximate  (EVI) j by 


(EVI) 


F ind  u^  e such  that 

a(uh>vh"uh)  - L(vh"uh)  VvheV 


Remark  1 . 7 : In  most  cases  it  will  be  necessary  to  approximate  also  a(*,*) 
and  L(*)  by  a^(»,*)  and  L^(»)  (usually  defined  - in  practical  cases  - from 
a(*,*)  and  !,(•)  by  a numerical  integration  procedure) . Since  there  is  nothing 
new  on  that  point  compared  to  the  classical  linear  case  we  shall  not  insist  more 
about  this  problem  for  which  we  refer  to  CIARLET  [33]. 

It  is  easily  proved  that 


Proposition  1 . 1 : Problem  (EVI).  has  a unique  solution. 
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1.5.3.  Convergence  results. 

Let  us  prove  now  the  following  convergence  theorem 


Theorem  1.3  : If  the  above  hypotheses  on  (V,  ).  . (K,  ),  hold  then 
n h h h 


(1.16)  lim  ||  u.  — u 1 1 = 0 

Vh-0 

where  u,u^  are  respectively  solutions  of  (EVI) ^ , (EVI) 


Proof  : (1)  Estimates  for  u, 
— r 

We  have  from  (EVI)., 


a(uh,uh)  < a(uh,vh)-L(vh-uh) 


which  implies  V v^  e that 


(1.17)  a|!uh||‘<  ||*||  IKHIIvJI  ♦ HlII,  ||»h||  + ||L||>||uh| 


From  (ii)  in  Sec.  1.5.1.  and  (1.17)  we  obtain 


ctlluj2  S||A||||uh||||rhv|M|L||J|rhv||+||L||J|uh||  VveX 


Now  take  vq  e X • Then  with  denoting  various  constants  depending  on  v but 
not  on  h,  we  have  from  (ii) 


r.v 
1 ho1 


< C 


V h 


which  implies 


(1.18)  a ||  uh]|  2 < C]  || uh||  + C2  Vh  . 

In  turn  (1.18)  implies  the  boundedness  of  (u,  ),  in  V. 

— h h 


(2)  Weak  convergence  of  (u^)^. 

We  can  extract  from  (u^)^  a subsequence,  still  denoted  by  (u^)^,  such  that 
(1.19)  u^  -*■  u weakly  in  V. 
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From  (i)  in  Sec.  1.5.1.  we  have 

(1.20)  u*  e K. 


From  (1.19)  and  property  (ii)  we  obtain  taking  the  limit  in 
a(uh,uh)  < a(uh,rhv)-L(rhv-uh)  VveX 

that 

(1.21)  lim  inf  a(u,,u,  ) < a(u  ,v)-L(v-u  ) V v e Y . 

u n h h n 

h-*0 

Then  we  observe  that 

alluh_u*i|2  - a(uh-u*,uh-u*)  = a(uh,uh)+a(u*,u*)-a(uh,u*)-a(u*,uh) 

implies  in  the  limit 

(1.22)  lim  inf  a(uh,u^)  > a(u  ,u  ) 

h -*•  0 

which  is  a (well-known)  weak  lower  semi  continuity  property. 
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We  have 


(3)  Strong  convergence  of  (u^)^. 

12 

Ot  II  uh~ u 1 1 < a(uh-u,uh-u)  = a(uh,uh)+a(u,u)-a(uh,u)-a(u,u1  ) - 

< a(uh>rhv)-L(rhv-uh)+a(u,u)-a(uh,u)-a(u,uh)  Vvey  • 

From  the  above  relation,  from  property  (ii)  and  from  the  weak  convergence 
of  (u^)^,  we  have  in  the  limit 

2 

a lim  sup  ||u^-u||  < a(u,v)-L(v-u)-a(u,u)  = a(u ,v-u) -L(v-u)  Vvcy  • 

But  since  X=K  we  also  have 

2 

(1.24)  a lim  sup  ||u,-u||  < a(u,v-u)-L(v-u)  VveK. 

Taking  v=u  in  (1.24)  it  follows  that 

2 

a lim  sup  ||u^-u||  ^0 

which  implies  that  lim  ||u, -u||  = 0 i.e.  the  strong  convergence. 

h+0 

Remark  1 .8  : Error  estimates  for  some  EVI's  of  the  first  type  have  been 
obtained  by  several  authors  (see  Sec.  2 for  more  details).  But  like  in 
many  non-linear  problems,  the  methods  used  to  obtain  these  estimates  are 
specific  to  the  particular  problem  under  consideration. 

This  remark  still  holds  for  the  approximation  of  EVI's  of  the  second  kind 
which  is  the  subject  of  the  next  sub-Section  1.6. 

1.6.  Internal  Approximation  of  (EVIjj. 

The  assumptions  on  V,  a(*,*),  L(*),  j(')  are  those  of  Sec.  1.2.  Furthermore 
we  assume  for  simplicity  that 

(1.25)  j(.)  is  continuous  over  V. 

I 

(In  Sec.  5 an  important  family  of  (EVI^s  for  which  j(»)  is  non  continuous 
will  be  considered  also) . 

-10- 
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1 .6.1.  Approximation  of  V . 

The  space  V is  approximated  by  the  family  (V  ) , being  a closed  subspace 
of  V(dim  c +=°  in  applications).  We  assume  that  (V^)  has  the  following 
property 

(i)  There  exist  If  cV,^=V  and  r,  : If  -*■  V,  such  that 
h h 

lim  ruv  = v strongly  in  V,  Vv  £ V . 
h-*0 


1.6.2.  Approximation  of  j(«). 

The  functional  j(*)  is  approximated  by  (j^)^  where 


(1.26) 


I is  convex,  1 . s .c. , uniformly  proper  in  h ; 

the  last  property  implies  the  existence  of  A : V •+  K,  linear 
and  of  )j  £ E such  that 

(1.27)  Jh(vh)  - Mv^)  + ^ VvheVh,Vh  ’ 

We  shall  assume  also  that  (j^)^  obeys 


and 


continuous 


(ii)  Lf  -*  v weakly  in  V then 

lim  inf  jh(vh)  £ j (v) 

(iii)  lim  j,  (r  v)  = j (v)  V v e V . 
h+0 

Remark  1 , 9 : In  all  the  applications  we  know,  if  j(’)  is  continuous,  then 
it  is  always  possible  to  construct  continuous  j^(*)  obeying  (ii),  (iii). 

Remark  1.10  : If  jh=j  yh,  then  (1.26),  (ii),  (iii),  are  automatically 
satisfied. 
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1.6.3.  Approximation  of  (EVI)^ 
We  approximate  (EVI)2  by 


(EVI)2h 


Find  u,  e V,  such  that 
n n 


1 a(uh*Vuh)+jh(vh)’ih(uh)  ~L(VV 

It  is  then  easily  proved  that 


V v e 
h 


Proposition  1.2  : If  the  above  hypotheses  on  (1/^)^*  (j^)^  hold  then  (EVI)^ 
has  a unique  solution. 

Remark  1.11  : Remark  1.7  of  Sec.  1.5.2.  still  holds  for  (EVI)  . 



1.6. A.  Convergence  results. 

Using  a variant  of  the  proof  of  Theorem  1.3  we  obtain  (see  GLOWINSKI  [53] 
for  more  details) 


Theorem  1 .4  : If  the  above  hypotheses  on  (V^)^,  (j^)^  hold  then 


lim  ||  u -u||  = 0 , lim  j (u,  ) = j(u) 
h-0  h h-0  h h 


where  u,uh  are  respectively  solutions  of  (EVI)2>  (EVI),^. 


2.  - SPECIFIC  EXAMPLES  AND  ERROR  ESTIMATES  FOR  CONFORMING  FINITE  ELEMENT 
METHODS . 

2.1.  Stationary  obstacle  problems. 

Problems  of  this  type  are  fairly  simple  but  provide  a good  mathematical  model 
for  several  important  applications.  Futhermore  obstacle  problems  are  those  for 
which  the  finite  element  approximation  error  analysis  is  the  most  achieved. 

2.1.1.  Formulation  of  a particular  obstacle  problem. 

N 

Let  ft  be  a bounded  domain  of  F with  a smooth  boundary  T = 3ft.  We  consider 
with  x = tx£^_|  and  ^ = {fix-  ^ the  Part^-Cu^ar  obstacle  problem 
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Find  ue  K such  that 


(2.1) 


Vu* V ( v-u) dx 


f (v-u) dx  Vv  e K , 


where  in  (2.1),  feL  ((2)  and  K is  defined  by 


(2.2)  K = { v t H ((2),  v > \p  a.e.  on  fi,  v|p  = g} 


with  4)  and  g given  functions  defined  respectively  on  Q and  T. 


2.1.2.  Physical  interpretation. 

2 

Assume  that  (2  c ; then  a classical  interpretation  of  (2.1),  (2.2)  is  that 
u represents  the  small  vertical  displacements  of  an  elastic  membrane  (2  under 
the  effects  of  a field  of  vertical  forces,  whose  intensity  is  given  by  f 
(f  represents  a surface  density  of  vertical  forces).  This  membrane  is  fixed 
on  its  boundary  T (u=g)  and  lies  over  an  obstacle , whose  height  is  given  by 
i{j(u  > i[i)  ; see  Fig.  2.1  for  a geometrical  description  of  the  phenomenon. 


Figure  2 . 1 

2.1.3.  Other  phenomena  related  to  obstacle  problems. 

Similar  EVI's  also  occur,  sometimes  with  other  type  of  boundary  conditions 
and/or  non  symmetric  bilinear  forms,  in  mathematical  models  for  the  following 

problems  : 

. I.ubrication  phenomena  (see,  e.g.,  CRYER  C 4 1 ] , [42],  MARZULLI  [67],  GLOWINSKI- 
LIONS-TREMOLIERES  [57,  ch.  2,  Sec.  5]  for  finite  difference  treatments  and  more 
references,  and  CAPRIZ  [27]  for  a discussion  of  the  mathematical  modelling). 


_ — , - q 

. Filtration  of  liquids  in  porous  media  (see  in  particular  BAL0CCH1  I J, 

L 2]  , [3],  C0MMINCI0L1  [35],  BAIOCCHI-BREZZI-COMMINCIOLI  L 5 ] , CRYER- 
FETTER  [43]  , BAI0CCH1-CAPEL0  [6  ] and  the  numerous  references  therein). 

. Tmu  dimensional  potential  flows  of  inviscid  fluids  (cf.  BREZIS-STAMPACCHIA 
[17]  , [18],  BREZIS  [11],  CIAVALDINI-POGU-TOURNEMINE  [34],  ROUX  [78]  and  the 
references  therein). 

. Wake  problems  (cf.  BOURGAT-DUVAUT  I"  1011  . 

This  list  is  far  from  complete  and  we  also  have  applications  in  Biomathematics, 
Economics,  Semi-conductors,  etc... 

2.1.4.  Interpretation  of  (2.1),  (2.2)  as  a free  boundary  problem . 

Let  define  from  the  solution  u of  (2.1), (2. 2) 


Sl+  = (x|xeSl,  u(x)  > (x)  } , 
(1°  = (x|xeS2,  u ( x)  = ip(x)}  , 


and  then 


Y - 9(1  n 9(1°  , 


u*  ‘ V 


Classically  (2.1),  (2.2)  has  been  formulated  as  the  problem  of  finding  v 
(the  free  boundary)  and  u such  that 


(2.3) 

(2.4) 

(2.5) 

(2.6) 


- Au  = f on  SI  , 
u = on  Sl°  , 
u = g on  T, 


U+K  = uo|y 


The  physical  interpretation  of  (2.3)-(2.6)  is  the  following  : (2.3)  means 
that  on  SI  the  membrane  is  strictly  over  the  obstacle  and  has  a purely  elastit 
behaviour  ; (2.4)  means  that  on  Sl°  the  membrane  is  in  contact  with  the  obstacle 
(2.6)  is  a transmission  relation  on  the  free  boundary. 
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r * 

In  fact  (2.3)— (2.6)  are  not  sufficient  to  characterize  u,  therefore  it  is 

necessary  to  add  other  transmission  properties  ; for  instance  if  t|i  is  smooth 

2 

enough  (let  say  \p  eH  (ft)),  we  should  require  the  "continuity"  of  Vu  on  y 
(we  may  ask  Vu  e H* (ft)  x H* (ft) ) . 

Remark  2 . I : This  kind  of  free  boundary  interpretation  holds  for  the  other 
examples  considered  in  the  sequel. 

2.1.5.  Existence,  Uniqueness,  Regularity  of  the  solution. 

Concerning  the  existence  and  uniqueness  of  a solution  of  (2.1), (2. 2)  we  can 
easily  prove 

1 I /2 

Theorem  2 ■ 1 : Assume  that  T is  smooth  and  that  ip  e H (ft) , g e H (T)  with 
ip|  p < g a.e.  on  T 

then  (2.1),  (2. 2)  has  a unique  solution. 

Remark  2.2  : The  above  theorem  holds  for  f e (h'((2))'  and  for  fairly 
discontinuous  tj>. 

Concerning  the  regularity  of  u let  us  recall  the  following  classical  results 
of  BREZIS-STAMPACCHIA  [ 1 9]  : 

If  T is  sufficiently  smooth,  if , for  p e ] 1 ,+°°  [ , f £ Lp(ft)  n (H* (ft)) ' , 
e W2  ’ P (ft)  , g = g|p  with  g e W2  ’ P (ft)  , then  u e W2 ’ p (ft) . 

Actually  the  above  results  have  been  refined  by  BREZIS  [12], [13]  and  very  sophis- 
ticated properties  of  the  solution  and  of  the  free  boundary  have  been  obtained 
by  Lewy-Stampacchia,  Brezis,  Kinderlehrer , Nirenberg,  Schaeffer,  etc... 

2.1.6.  Finite  element  Approximations  of  (2.1), (2. 2).  (I)  Piecewise  linear 
Approximations . 

We  consider  in  this  subsection  conforming  finite  element  approximations  of 
order  one  of  the  obstacle  problem  (2.1), (2. 2).  Piecewise  quadratic  approximations 
are  considered  in  Sec.  2.1.7  and  non  conforming  approximations  of  mixed  type  in 
Sec.  4. 
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We  assume  for  simplicity  that  £1  is  a bounded  polygonal  domain  of  F . 

We  assume  also  that  ijj  e H1  (£2)  n C°(£2)  , g e H 1 (D  n C° (D  . We  introduce  then 
a standard  triangulation  of  £7  such  that 

U T = £2  , 

T£<Gh 

with  as  usual,  h = length  of  the  largest  side  of  C!  . 


Let  define  now 


E,  = {Pe  £2,  P vertex  of  T e } 
n n 


E,  = { P e E,  , P^D  = E,n  £2 . 
h h h 


We  approximate  then  H (£2)  and  K by  respectively 


(2.7)  Vh  = (vh€  C°(S),  vh|Te  P,  V T e 'C^  , 

(2.8)  = (vhe  Vh,  vh(P)  > 4»(P)  VPe  Zht  vh(P)  = g(P)  VPe  Ehn  D 

where  in  (2.7)  and  the  sequel  (for  ks  0)  P^  = space  of  polynomials  in  two 
variables  of  degree  S k. 


Finally  we  approximate  (2.1),  (2.2)  by  (2.8)  and 


Find  u^e  such  that 


(2.9) 


* , 

,Q  Vuh*V(vh-uh)dx  > Jnf(vh-uh)dx  Whe  1^. 


Proposition  2 . 1 : The  approximate  obstacle  problem  (2.8),  (2.9)  has  a unique 
solution. 

Concerning  the  convergence  of  the  approximate  solutions  as  h + 0,  we  refer 

to  GLOWINSKI  [53,  Chapter  4,  Sec.  2]  for  the  case  where  u is  not  very  smooth. 

Below  we  shall  briefly  consider  the  derivation  of  error  estimates,  in  the 
1 . ~ 2 

H (£2)  norm,  if  u,i|/,  g cH  (£2).  We  follow  very  closely  the  analysis  of  BREZZI- 


HAGER-RAVIART  [21,  Sec.  4]  (in  which  a( 


u,v)  = 


Vu*Vv  dx  + uv  dx)  . 
£2  J£2 
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To  obr.tir.  these  error  estimates  we  need  the  following 


Lemma  2.1  : Assume  that  ueH  (SI)  ; then 


(2.10)  - Au-f  -0  a.e.  on  SI, 


(2.11)  (-Au-f ) (u-tjj)  = 0 a.e.  on  SI. 


Proof  : see  BREZ1S  [14 ]. 


Lemma  2 , 2 : Let  u and  be  respectively  solutions  of  (2.1),  (2.2)  and  (2.8)  , 
(2.9).  Then 


(2.12)  a(u  -u,u,-u)  $a(u,'u,v,-u)+a(u,v,-u,)- 

h h h h h h 


f(vh-uh)dx  VvheKh. 


Proof  : Following  BREZZI-RAGER-RAVIART  [21  , Theorem  2.1]  we  have  Vv^ e 


a(uh~U ’uh_U)  = a(uh~u»vh"u)+a(uh~u,uh_Vh^ 


\ = a(uh-u,vh-u)+a(u,vh-uh)-  Jftf(vh-uh)- 


.f(vh~uh)dx'a(uh*vh~uh) 


Since  obeys ( 2 . 9),  we  have 


f(Vuh)dx~a(uh-vh'uh)  -°  Vvh£Kh 


which  combined  with  the  above  equation  implies  (2.12). 


We  prove  now 

Theorem  2.2  : f £ I." (SI)  ,ip  £ H^(S1)  , g = g 1 p . g £ H^(S1)  and  if  the  angles  of 

are  bounded  below  by  8^  > 0 , independent  of  h , then 

(2.13)  !|u,-u||  = 0(h)  , 

n H (SI) 

where  u and  u^  are  respectively  solutions  of  (2.1),  (2.2)  and  (2.8),  (2.9). 
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Proof  : We  follow  again  BREZZI-HAGER-RAV1ART,  loc.cit.,  Theorem  4.1  (see  als 
FALK  [48]).  We  have  from  Green's  formula 


(2.14)  a(u,v)  = Vu*Vv  dx  = - Auvdx+  4~“  v dT  V v e h'  (12)  . 


Since  v,  -u,  e H ((2)  Vv,  e K,  it  follows  then  from  (2.14)  that 
h h o h h 


(2.15)  a(u,vh~uh)-  J f(vh~uh)dx  = j (-Au-f) (vh~uh)dx  V vh € K^. 


Let  be  the  operator  of  V^-interpolation  on  E^.  Then  since  SlcR  we  have 
H“((2)  c C°(f2)  and  ueH2((2),  u > i|)  on  (2  imply 


(2. 16)  e K^. 


Taking  v^  = tt^u  in  (2  . 1 2)  , (2 . 1 5)  we  obtain 


(2.17)  a(u^-u,u^-u)  < a(u^-u,TT^u-u)+J  (-Au-f ) (tt^u-u^) dx 


Observe  that 


(2.18)  Trhu-uh  = (7rhu~u)  + ^_Trh^+(u“'lJ)  + ('rThl/J-uh)  • 


Let  w = -Au-f  ; then  weL  ((2)  and  from  (2.18) 


(2.19) 


w(tt,u-u,  )dx  = w(TT,u-u)dx  + w(i/j-n,ip)  dx  + w(u-’^)dx  + 

(2  JS2  J(2  '(2 


+ Jfiw(V~uh)dx  • 


Since  Lemma  2.1  holds  we  have  w^O  a.e.  and  w(u-i|i)  = 0 a.e.  ; moreover  since 
u^eK^  we  have  TT^lJi-u^  S 0 on  (2.  It  follows  then  from  (2.19)  that 


(2.20)  w(tt  u-u  )dx  <||w||  „ ( ||tt  u-u||  . + ||tl<HHI  2 ) 

h h h IT  ((2)  h LZ((2)  h L^(Q);' 


Since  u.iJj  e H ((2)  we  have  (since  the  angle  condition  holds) 


(2.2!)  ||  tt  u-u  1 1 < Ch  ||u||  . ,||  7T  u-u  |l  < Ch2  ||u|| 

H (£2)  H (52)  Lz((2)  H ((2) 


i 


(2.22) 


< C h 


H ((2) 


H2  (SI) 


L2((2) 


Ch2  ||  ^1! 


H2  (Ti) 


where  the  C’s  are  independent  of  h,u,i]/  . 

It  follows  then  from  (2 . I 7) , ( 2 . 20) -( 2 . 22)  that 

<2.23)  I V“l  , ,0  - DO) 
where 

I Vl  1 ({n  | Vv|  2dx) 

(and  where  below  ||v|j  = )|v||  )• 

H (52) 

To  estimate  ||  u^— u | ] j ^ we  observe  that,  since  (2  is  bounded 
(2.24)  llv||]  ft“^lvli  f2  ^ v 6 H^((2)  , C independent  of  v . 
It  follows  then  from  (2.24)  and  from  u,  —it.  u e H*  (52)  that 


|llvulll,52  ^Huh_7Thull  t a+  l|7rhu"u|1,  n - 

sCl  Wli.n + H V-"!!,  n - 

“ Cl  uh_ul  1 ,52  + Clu"7rhull,52+  Ku'uHl,52  * 
Then  (2:13)  follows  clearly  from  (2 .2 ! ) , (2 . 23)  , (2 . 25)  . 


To  our  knowledge  the  first  0(h)  estimates  for  ||u^-u||  j have  been  obtained, 

for  piecewise  linear  approximations,  by  FALK  [46]  and  then  M0SC0-STRANG  [72]  . 

These  works  have  been  followed  by  FALK  [47],  [48]  (see  also  GLOWINSKI-LIONS- 

TREMOLIERES  [57,  Ch.  1],  CIARLET  [33]).  In  our  opinion  one  finds  in  FALK 

[48]  the  most  complete  analysis  for  piecewise  linear  approximations,  since 

it  also  considers  non  convex  and/or  non  polygonal  fl.  The  problem  of  obtaining. 

2 

via  a generalization  of  the  Aubin-Ni tsche  trick,  L -estimates  (of  optimal 
order)  is  not  completely  solved  yet  ; however  for  some  partial  results  in  that 
direction  see  NATTERER  [76],  M0SC0  [70]  and  the  references  therein. 
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To  conclude  with  piecewise  linear  approximation,  let  us  mention  that  under 

suitable  hypothesis  BAIOCCHI  [ 4l  (resp.  NITSCHE  [77])have  obtained  for  the 

obstacle  problem  ||u  -u]|  = CKh"  £)  , e>0,  arbitrarily  small  (resp. 

" L (ft) 

lluh'u|>  =0  = 0(hZ|Log  h|)). 

L (ft) 

2.1.7.  Finite  Element  Approximations  of  (2.1),  (2.2).  (II)  Piecewise  quadratic 
approximations . 

O 

With  £.  , E,  as  in  Sec.  2.1.8  define 
h h 

E^  = (Peft  , P midpoint  of  a side  of  T e , 


E^  = {PeE^  , P i T)  , 


O o o 


E"  = , E"  = E,  u . 

h h h h h h 


We  approximate  H (ft)  and  K by 

(2.26)  vh  = {vh  eC°(ST)  , vhlTeP2  VTe^}  , 

(2.27)  j = (vhe  Vh  , vh(P)  a i|)(P)  V P e vh(P)  = g(P)  VPeZjJ  n T}  , 

(2 . 27)  2 kJ  = (vheVh  , vh(P)  2l|»(P)  VPeE^  , vh(P)  = g(P)  VPeE^nT}  . 

2 

We  observe  that  in  the  condition  v^(P)  - ip(P)  is  only  required  on  the  side 
midpoints . 

We  approximate  the  obstacle  problem  (2.1), (2. 2)  by 


F ind  u^  e such  that 


(2.28) . 


( L VVV(vh~Uh)dx  " Lf(vh“Uh)dx  VvheKi  * 


where  i=l ,2. 

Proposition  2.2  : The  approximate  problems  (2.28)^  have  a unique  solution. 

Concerning  the  convergence  of  the  approximate  solution  we  refer  to  GLOWINSKI 

[53]  where  lim  ||u^-u||.  = 0 is  proved,  for  i=l,2,  assuming  the  usual  angle 

h-»-0  h >“ 

condition.  The  error  estimates  analysis  is  much  more  complicated  that  with 
k=l , and  we  refer  to  BREZZI-HAGER-RAVIART  [21 3 where  under  suitable  assump- 


t 


t ions  on  f9u9\p9g  and  the  free  boundary,  one  proves  that  Vi  = l,2 


i uh  ull  | ~ 0(h  b)  , e>0  arbitrary  small, 

H1  (ft) 

if  the  angle  condition  holds. 

2* l *®*  Concluding  Remarks.  Further  comments. 

Non  linear  obstacle  problems  have  been  considered  by  WH1TEMAN-N00R  [83]. 
Concerning  the  numerical  solution  of  the  approximate  obstacle  problems  we 
refer  to  Sec.  3.1  where  several  methods  of  solutions  are  described. 


2 • 2 . The  elasto-plastic  torsion  problem. 

The  problem  we  shall  consider  in  this  section  is  more  complicated  than  the 
obstacle  problem  of  Sec.  2.1.  It  is  related  to  the  elasto-plastic  torsion 
of  a cylindrical  bar  of  infinite  length.  It  is  a fairly  simple  plasticity 
problem  but  it  is  of  great  interest  from  both  theoretical  and  numerical 
points  of  view. 

2.2.1.  Formulation  of  the  continuous  problem. 

2 

Let  ft  be  a bounded  domain  of  R with  a smooth  boundary  f = 3ft  . We  consider 


Find  ueK  such  that 


(2.29) 


Vu*V(v-u)dx  > 


f (v-u)dx  V v e K, 


ft 


where  f e L (ft)  and 


(2.30)  K = {v  e IT  (ft)  , | Vv|  <1  a.e.}  . 


We  recall  that 


Hq (ft)  = JH&)  H (j,  ) = { v | v e H ' (ft)  , v-0  on  8ft}, 


Remark  2.2  ; Since  the  bilinear  form  a(*,*)  occuring  in  (2.29)  is  symmetric 
fa(u,v)  = j 7u»Vv  dx) , (2.29)  is  indeed  equivalent  to  the  minimization  problem 
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Find  u t K such  that 


(2.31) 

[ J(u)  < J ( v ) V v e K, 


where  J(v) 


|Vv|dx 

ft 


f v dx . 

ft 


2.2.2.  Physical  interpretation. 

Let  consider  an  infinitely  long  cylindrical  bar  of  cross  section  ft  , where 
is  simply  connected.  Assume  this  bar  is  made  of  an  isotropic,  elastic, 
perfectly  plastic  material  whose  plasticity  yield  is  given  by  the  Von  Mises 
criterion  (see  DUVAUT-LIONS  [45,  Ch.  5]  , for  a general  discussion  of  plasti- 
city phenomena).  Starting  from  a zero-stress  initial  state,  an  increasing 
torsion  moment  is  applied  to  the  bar.  The  torsion  is  characterized  by  its 
torsion  angle  per  unit  length  C.  It  follows  then  from  the  Haar-Karman  prin- 
ciple that  the  stress  field  can  be  obtained  through  the  solution  of  the 
following  variational  problem 


(2.32) 

which  is 


Min  { y | Vv | ^dx  - C 
veK  ft 

a particular  case  of 


v dx} 

Jq 

(2.29)  , 


(2.31), 


with  f=C. 


The  stress  vector  0 in  ft  is  obtained  from  u by  0 = Vu.  Hence  u appears  as 
a stress  potential. 

Remark  2.3  : If  ft  is  not  simply  connected,  the  formulation  of  the  elasto-plast . 
problem  has  to  be  modified  and  we  refer  to  GLOWINSKI-I.ANCHON  [56],  GLOWiNSKl 
[53,  Ch.  4,  Sec.  3.2]  for  the  new  formulation. 


2.2.3.  Existence  and  Uniqueness  results.  Regularity  and  further  properties. 

The  condition  of  Sec.  1 being  fulfilled  we  can  apply  Theorem  1.1.,  then 

Proposition  2.3  : Problems  (2.29),  (2.31)  have  a unique  solution. 

For  the  proof  we  refer  to  GLOWINSKI-LIONS-TREMOLIERES  [57,  Ch . 31  . 

Concerning  regularity  properties  it  follows  from  BREZIS-STAMPACCHIA  Ip; 
that  if  3ft  is  smooth  (or  ft  convex)  and  if  feL^(ft)  with  p 2 then  the  solution 
u of  (2.29),  (2.31)  satisfies 
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u e W2,P(ft)  n K. 


If  in  particular,  f = const,  (as  in  (2.32))  then  u e_  W *^(ft)  for  p arbitrary 
large.  If  for  example  £2  is  a disk  and  f - const,  then  for  f large  enough 
u e w2  ’“(ft)  (1HS(!1)  V s < | , but  u 4 C 2 (ft)  , u 4 HJ  (ft)  . 

Remark  2.4  : If  f = const.,  BREZIS-SIBONY  [20]  have  proved  that  the  solution 
u of  (2.29),  (2.31)  is  also  the  unique  solution  of  the  two-obstacles  problem 


(2.33)  Min  i |Vv|2dx  - C v dx) 
veK*  Jft  'ft 


where 


K = (ve  HQ(ft)  , |v(x)|  < 6(x,D  a.e.} 


with  6(x,D  = distance  of  x to  F = 3ft. 

Remark  2.5  : For  the  free  boundary  asnprf  of  the  elasto-plastic  problem  we 
refer  to  e.g.  GLOWINSKI-LIONS-TREMOLIERES  [57].  Actually  SHAW  [80]  has  nume- 
rically solved  (2.32)  as  a free  boundary  problem  using  finite  difference 
approximations . 

2.2.4.  Finite  element  approximations  of  (2.29),  (2.31). 

We  assume  ft  polygonal , then  we  define  as  in  Sec.  2.1.6  and  we  introduce 


5„t,  - NkeC°<a)  , vjr  - 0 , vh|TeP|  VTeCh)  , 


5,  ' KeVh  • lWtJ  5 1 * '""'oh  ■ 


Then  we  approximate  (2.29)  by 


Find  u^eK^  such  that 


(2.34) 


Vu,  •V(v, -u,  )dx  5 f(v,-u,)dx  V v.e  K.  . 

Jo  h h h'  Jo  h h h ^ 


It  is  clear  that  (2.34)  has  a unique  solution. 


Remark  2.6  : Since  Vv^  is  piecewise  constant  , the  condition  v^  e K^  amounts 
to  Card  ('£’  ) quadratic  constrains  (|Vv^|^5  1 on  T.VTeC^).  If  instead  of 
a piecewise  linear  approximation,  one  uses  a piecewise  quadratic,  requiring 
Vv^|  < 1 a.e.  then  would  amount  to  3 CardCJ^)  quadratic  constraints 

(see  GLOWINSKI-LIONS-TREMOLIERES  [57,  Ch.  3]  for  more  details). 

Remark  2.7  : The  numerical  analysis  of  (2.29)  via  (2.33)  is  done,  if  f = const, 
in  GLOWINSKI-LIONS-TREMOLIERES,  loc.  cit.,  Ch.  3. 

H*  (ftl  , co  _ 

Convergence  analysis  : Since  .#(ft)  nK  ov“  = K (where  ,£(ft)  = {v  e C (ft), 
v has  a compact  support  in  ft})  we  can  prove  using  the  general  approximation 
results  of  Sec.  1 that 

(2.35)  lim  ||u  -u||  = 0 

h-H)  ’ 

if  the  angle  condition  holds.  For  the  proof  of  (2.35)  and  of  the  above  density 
result  see  e.g.  GLOWINSKI-LIONS-TREMOLIERES,  loc.  cit.,  Ch.  3 and  GLOWINSKI 
[53,  ch.  4,  Sec.  3]. 

Moreover  if  f c L^(ft)  and  ueW2’P(ft)  with  p>2,  it  is  proved  in  FALK  [463 
(see  also  GLOWINSKI  63,  Chap.  4,  Sec.  3]  ) that 

(2.36)  II  uh-u  1 1 1 ,ft  = 0(h,/2_1/P)  . 

In  FALK  [46l  one  also  considers  the  case  where  ft  is  non  polygonal. 

2 

Remark  2.8  : If  ft  c R then  f eL  (ft)  implies  that  ||u^-u||j  ^ = 0(h),  instead 
of  (2. 36). This  result  is  related  to  the  fact  that  in  the  monodimensional  case 
the  piecewise  linear  interpolate  of  v e K is  still  in  K,  which  is  no  longer 
true  in  R . 

The  iterative  solution  of  the  approximate  problems  is  discussed  in  Sec.  3.2. 

2.3.  Flow  of  a Bingham  medium  in  a pipe. 

We  have  considered  in  the  two  previous  sections,  examples  of  EVI's  of  the  first 
kind.  In  the  present  section  we  shall  discuss  an  EVI  of  the  second  kind  related 
to  the  flow  of  a Bingham's  viscous-plastic  medium  in  a pipe.  This  section  follows 
GLOWINSKI  [54 ] and  GLOWINSKI-LIONS-TREMOLIERES  [57,  Ch.  5]  (see  also  DUVAUT- 
LIONS  [45,  Ch.  6]  for  a more  precise  mechanical  interpretation). 
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2.3.1.  Fo rniulation  of  the  continuous  problem. 

2 

Let  [1  be  a bounded  domain  of  R“  with  a smooth  boundary  T = Let  define 

j(')  by 

(2.37)  j(v)  = |Vv|dx  ; 

h 

j(*)  is  Lipschitz  continuous  but  not  differentiable. 

2 

Let  consider  now  the  following  EVI  of  the  second  kind  (with  feL  (12))  : 

!F i nd  u e (12)  such  that 

y Vu*V(v-u)dx  + gi  (v)-gj  (u)  > f(v-u)dx  VveH*(12)  , 

h J12  0 

which  is  equivalent  to 
— 

(2.39)  Min  {y  | Vv | dx  + gj (v)  - fv  dx}  . 

veH^ (Q)  hi  h 

Assuming  that  y >0  and  g 5 0 if  follows  from  Theorem  1.2  of  Sec.  1.4  that 

Proposition  2.4  : The  two  equivalent  problems  (2 . 38) , (2 . 39)  have  a unique 

solution . 

2.3.2.  Mechanical  Interpretation. 

If  f = const.  = C ( C>0  for  example)  it  follows  from  LIONS-DUVAUT  [ 45  , 

Ch.  6]  that  (2 . 38) , (2 . 39)  models  the  laminar  stationary  flow  of  a Bingham's 
viscous  plastic  fluid  in  a cylindrical  pipe  of  cross  section  fi,  with  u(x) 
the  velocity  at  x.  The  above  constant  C is  the  linear  decay  of  pressure 
and  y,g  are  respectively  the  viscosity  and  the  plasticity  yield  of  the 
medium . 

The  above  medium  behaves  like  a viscous  fluid  (of  viscosity  y)  in 
f2+  = {xjxefi  , |Vu(x)  | >0} 


and  like  a rigid  medium  in 


Find  u,  c V , such  that  V v,  £ V , 
h oh  h oh 


(2.43) 


y VVV(vh-uh)dx  + g |Vvh|dx  - g |Vuh|dx  5 f(vh-uh)dx  . 


The  approximate  problem  (2.43)  has  clearly  a unique  solution.  Concerning 
the  convergence  of  u^  to  u as  h •*  0 we  have 

Theorem  2.3  : Assume  that  the  angles  of  are  uniformly  bounded  from 

below  by  0 >0,  as  h -*■  0,  then 

— o 


(2.44)  lira  ||  u h-u  ||  , o=  °* 


If  furthermore  u cH  ($7),  then 


(2.45)  II uh-u ||  , ^ = 0(h  ) . 

Proof  : We  follow  GLOWINSKI  [54].  Taking  v^=0  in  (2.43)  we  obtain 

(2-46)  kl?  o-  illfH  2 Nuhll  2 Vh  • 

f 2 1/2 

Since  is  bounded  |v|,  n = ( | |Vv|  dx)  defines  on  H (ft)  a norm  equi- 

J q 0 

valent  to  ||v||j  We  have  moreover 

1IMI  2 5 ~TT2  Mi  n Vve  » 

(2.47)  V ’ 

[ Xq  = smallest  eigenvalue  of  -A  over  H^(ft). 


It  follows  then  from  (2 . 46) , ( 2 . 47)  that 


(2-48)  Kli  ft-  — T77  HfH  2 

yX  LZ(ft) 

o 

In  other  respect  we  do  have 


% Vuh*V(vh'Uh)dx  + 8 " jnf(vh-Uh)dx  V V Voh  • 

w|  Vu*V (u^-u)dx  + gj(uh)-gj(u)  > f (uh~u) dx  , 


and  hence  by  addition  we  obtain 


(2.49)  p|uh-u|“  ^<gj(vh)-gj(u)+yj  Vuh*V(vh-u)dx-  j f (vh~u) dx  V vh e Vh . 


From  (2.48) ,(2.49)  and 


(2.50)  j (v) 


|Vv|dx  < 
■ o 


meas . (Q.)  |v|  0 VveH  (ft) 

* f o l O 


we  obtain 


(2-51)  luh"ul  i - y (g  / meas-  (n)  + II  f HL2(n))  I Vh_U  I 1 ,ft  Vvh  6 Vh 

Let  <$>  e i)(ft)  ; we  denote  by  TT^tJ)  the  V^- interpolate  of  <J>  on  i.e. 


VeVoh  ’ 


tt^4>  (P)  = 4> ( P)  V P vertex  of 


Since  the  angle  condition  holds  we  have 


(2.52)  | it .<M>|  . < C||<fr||  . h V4>  e MW 

n i,u  h (ft) 

with  C independent  of  h and  <{>. 


From  (2.51) ,(2.52)  and  from  the  triangular  inequality  we  obtain,  taking 
vh  = TT^if)  in  (2.51) 


(2.53) 


IV^I.ft  " y (g 


V 4>  € ,£(ft) . 


+ 7 ==  llfll  2 )(|4'-u|,  o+CH(J,ll  2 h) 

'f  Xo  L^(ft)  ’ IT  (a) 


Since  £(0.)  is  dense  in  H_(i2),  (2.53)  implies  clearly  (2.44).  To  prove  (2.45) 

® 2 1 o 

we  use  directly  (2.51)  ; indeed  if  u e H (ft)  n 1)^(0)  , then  ueC  (ft)  and  tt^u 

can  be  defined.  We  have  furthermore 


(2.54)  Ku-u|  . Q<  C h ||  u 1 1 

’ H (ft) 

Then  taking  v^  = tt^u  in  (2.51)  it  is  clear  that  (2.54)  implies  (2.45).  Q.E.D. 
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Remark  2.9  : Assuming  some  extra-regularity  assumptions  (satisfied  if  12 
is  a disk  and  if  f = const.),  it  is  shown  in  GLOWINSKI  [53,  Ch.  5],  [ 54] 


Huh“ull  i ,n=  0(1 108  hl  h)  • 

Iterative  methods  for  solving  (2.38)  and  (2. A3)  are  described  in  Sec.  3.3. 


2.4.  Error  estimates  of  optimal  order  for  the  elasto-plastic  torsion  and 
Bingham's  flow  problems  via  a new  formulation. 

In  this  Section  we  shall  briefly  describe  some  of  the  results  of  FALK-MERCIER 
[49]  , who  via  a new  variational  formulation  of  the  problems  in  Sec.  2.2, 

2.3  have  obtained  error  estimates  of  optimal  order. 

2.4.1.  The  elasto-plastic  torsion  problem. 

We  consider  again  the  elasto-plastic  problem  of  Sec.  2.2.  We  assume  that 
12  is  simply  connected  (which  corresponds  precisely  to  the  physical  problem) . 


Proposition  2,5.  : Assume  that  12  is  simply  connected,  then  the  variational 
problem  (2.29)  is  equivalent  to 


Find  p £ A n H such  that 


(2.55) 


p’(q-p)dx  S <j)*(q-p)dx  Vq  e AnH 

12  J12 


3<t>j  9 4*2 

where  <|>  = (<j>j  ,<[),,}  is  any  solution  of  f = and  where 

(2.56)  A = (q  £ L2 (12)  *L2(12)  , |q|  < 1 a.e.  on  12} 

(2.57)  H = {q  £ L2(12)  xl2(12)  , [ q*Vw  dx  = 0 Vw£  H ' (12) } . 


The  solutions  u of  (2.29)  and  p of  (2.55)  are  related  b) 


<2-58>  p = • 

2 I 


4 


Remark  2.10  : It  is  clear  that  q e H is  equivalent  to 


l q e L2(P.)  x L2(p)  , 

(V*q  = 0 a.e.  on  P,  q*n  = 0 a.e.  on  9P 

where  n is  the  outward  unit  normal  vector  at  P. 


Remark  2.11  : If  £ = const.  = C,  one  can  take 


<j>  | = Cx2  , <|>2  = 0 


If  f is  not  constant  and  9P  smooth  or  P convex  it  is  always  possible  to 

. ~ "TT^s  , . , 3<fr|  3<t> 2 


construct  from  f £l  (P)  , $|,<J>2  £H  (P)  such  that  f = 


9x. 


9x , 


For  instance 


we  solve 


Au  = f on  P, 
o 


u tH  (P) 
o o 


2 1 

which  produces  u eH  (P)  OH  . It  suffices  then  to  take 
o o 


<p,  - 


8u 

3^  ’*2 


8u 
c 

8x, 


The  approximate  problem  : For  simplicity  we  assume  that  Q is  a polygonal 

bounded  convex  domain  of  ; let  tr.  be  a triangulation  of  P like  in  the 

" . I ” 2 2 

above  sections.  We  approximate  H (P) , L (P)  xL  (P)  and  H by  respectively 


Vh  = {vhec°(P)  ,vh|TeP|  VTety, 

Lh  = {qh  €L2(P)  xL2(P)  , qh|T£F2  VTe^}  , 


Hh  = {qh£V 


qh*^Wh  dx  = 0 V wu  tVL)  . 


h hJ 


Then  we  approximate  (2.55)  by 


Find  p,  e A (i  H,  such  that 
r h h 


(2.59) 


ph*(qh_ph)dx  2 


i‘(qh"Ph)dx  V qh  € A n Hh‘ 
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The  approximate  problem  (2.59)  has  a unique  solution  and  it  is  proved  in 
FALK-MERCIER  [49]  that  if  some  convenient  properties  of  u and  p hold  together 
with  the  usual  angle  condition,  then 


(2.60) 


lPh"pi 


L2(ft)x  L2(Q) 


0(h)  . 


If  for  example  f = const.,  which  correspond  to  the  physical  problem,  then 

(2.60)  holds.  We  refer  to  FALK-MERCIER,  loc.  cit.,  for  more  details. 


Concluding  remark  : In  conclusion  by  a change  of  formulation  an  error  estimate 
of  optimal  order  have  been  obtained.  It  seems  however  that  (2.59)  is  more 
complicated  to  solve  numerically  than  (2.34). 


2.4.2.  The  Bingham's  flow  problem. 

We  consider  now  the  Bingham’s  flow  problem  of  Sec.  2.3.  Then  we  have 
Proposition  2,6  : The  variational  problem  (2.38)  is  equivalent  to 


Find  p e H such  that 


(2.61) 


p*(q-p)dx  + gj(q)  - gj(q)  > 4>*(q-p)dx  VqeH 

n Jo 


where  <p  and  H are  like  in  the  statement  of  Prop.  2.5  and  where 


j(q)  = 


n 


I q I dx  . 


The  solutions u of  (2.38)  and  p of_  (2.61)  are  related  by 


/ 9u  _ 9u  i 

P = ldx2  ’ " 3x,  ' ' 


The  approximate  problem  : With  Vh,Lh,Hh  as  in  Sec.  2.4.1.  we  approximate 

(2.61)  by 


Find  p,  e H,  such  that 
rh  h 


(2.62) 


fQ  ph*(Vph)dx  + 8j(qh)  ' 8J(Ph)  2{^*(qh-Ph)dx  Vqh€Hh- 
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The  approximate  problem  (2.62)  has  a unique  solution  and  under  suitable 
assumptions  on  u and  p it  follows  from  FALK-MERCIER  [A9]  that 

l|Ph"Pll  2 2 0(h) 

" l/(ft)xl/(ft) 

(we  still  assume  that  the  angle  condition  holds) . 

The  concluding  remark  of  Sec.  2.4.1.  still  notds  for  (2.62). 

2.5.  Further  problems. 

We  have  not  considered  in  this  paper  the  numerical  analysis  via  finite 
elements  of  problems  like 


Find  u e K such  that 


(2.63) 


with 


Vu*V(v-u)dx  > 


f (v-u)dx  Vv  e K 


ft 


K = {veH  (ft)  ,v|  = g , v 2 ip  a.e.  on  T.) 

r 


where  T and  T.  are  such  that  T n T,  = 0 , F u T.  = 9ft,  or  like 

n i o 1 ~ ’ 


W/  X . 

o 1 


(2.64) 


with 


Find  u e H (ft)  such  that 


Vu»V(v-u)dx  + j (v)  - j(u)  2 


f(v-u)dx  VveH  (ft)  , 


ft 


j(v)  = g 


v dr 


9ft 


For  finite  element  approximations  of  these  problems  by  conforming  methods 
we  refer  to  GL0W1NSKI-LT0NS-TREM0LIERES  [57,  Ch.  4]  , GLOWINSKI  [53,  Ch.  4], 
SCARPINI-VIVALDI  [79],  MOSCO  [70],  BREZZI-HAGER-RAVIART  [21],  etc... 
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- ! IT  RAT  1 VI.  METHODS  FOR  SOLVING  THE  APPROXIMATE  VARIATIONAL  INEQUALITIES 
OF  SEC . 2. 

In  this  section  we  shall  give  some  brief  indications  on  the  actual  solution 
of  the  approximate  problems  of  Sec.  2 by  iterative  methods.  For  more  details 
and  ether  applications,  we  refer  our  readers  to  GLOWINSKI-LIONS-TREMOLIERES 
I 571  , TREM0L1ERES  181],  GLOWINSKI  [55],  MERCIER  [69],  the  references 
therein  and  the  references  below.  In  fact  the  methods  to  follows  are  closely 
related  to  Non  Linear  Programming. 

3.1.  Iterative  solution  of  the  obstacle  problem. 

3.1.1.  Orientation 

The  obstacle  problem  of  Sec.  3.1  and  its  variants  (possibly  involving  non 
symmetric  a (•>•))  may  be  solved  by  several  methods.  We  shall  concentrate 
on  S.O.R.  with  truncation  (see  Sec.  3.1.2),  duality  (see  Sec.  3.1.3)  and 
give  some  indication  on  the  use  of  penalty  methods  in  Sec.  3.1.4. 

Actually  solution  of  discrete  variational  inequalities  based  on  the  so-called 
Complementarity  Methods  have  been  studied  these  last  years  ; in  that  direction 
we  refer  to  e.g.  COTTLE  [38],  [39],  COTTLE-GOLUB-SACHER  [40],  MOSCO-SCARPINI 
i 71  i . In  our  opinion  these  methods  are  less  effective  (at  least  for  most 
discrete  EVI's)  than  the  methods  of  Sec.  3.1.2,  3.1.3  and  will  not  be 
considered  here. 

3.1.2.  Methods  of  S.O.R.  with  truncation. 

These  methods  have  been  by  far  the  most  popular  for  solving  discrete  obstacle 
problems  and  are  in  our  opinion  the  simplest  to  program  and  the  most  economi- 
cal in  term  of  computer  storage. 

The  various  discrete  obstacle  problems  we  have  discussed  in  Sec.  2.2  are  in 

fact  particular  cases  of 

(3.1)  Min  1 y (Ay ,y)-(b,y) } 
yrC 

N 

where  in  (3.1)  ,^y  = {yj  , . . .y^}  eH  , A is  a NxN  positive  definite  symmetric 

matrix,  fx,y)  = 7 x.y.  , bcK^  and 
i=l  1 1 ~ 
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(3.2)  C = { y e F‘  , a . < y . < b . V i = 1 , . . . N } 

i iii 


with  i b.  Vi=l,...N.  Some  of  the  a.'s  (resp.  b.'s)  can  possibly  be  equal 

to  - 00  (resp  . + 00  ) . 

If  A = (ajj)  j ^ j < then  a typical  S.O.R.  + truncation  algori  ttrni  is 

(3.3)  x0£R^  , arbitrary  given, 

n . n+1  - n 

assuming  x known,  one  computes  x from  x , component  by  component  as 
follows 


For  i=l,...  N,  compute 


(3. 

.4). 

n+1/2  _ 1 

i cL  . 

(b 1 a.  -xn+ 1 - 7 a . . 

ii 

1 j<i  1J  3 j>i  1J 

(3. 

,5)i 

n+l  _ , n 

x . = P . (x . 

1 11 

. n+1/2  n.. 

+ w(xi  -x.)) 

with  Pi(yi)  = sup  (a^ , inffb-.y^). 


About  the  convergence  we  have 


Theorem  3 . 1 : If  0<U)<2,  then  Vx°eIlN, 


lim  x = x 

n-M-  00  ~ 

where  x is  the  unique  solution  of  (3.1). 

For  the  proof  of  Theorem  3.1  we  refer  to  CRYER  [42],  CEA-GLOWINSKI  129], 
COMINCIOLI  [36],  GLOWINSKI-LIONS-TREMOLIERES  [57,  Ch.  2].  In  some  of  the 
above  references  generalization  to  block  algorithms  in  Hilbert  spaces  are 
also  considered. 


Remarks  on  the  choice  of  to  : It  appears  that  the  optimal  value  w of  uj 
r r opt 

(i.e.  this  giving  the  fastest  convergence  for  a given  norm)  is  a function 
of  C and  b.  Therefore  for  the  discrete  obstacle  problems  discussed  above 
it  will  be  a function  of  f,g,^. 
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In  practice  several  strategies  can  be  used  ; one  can  either  use  the  optimal 
u)  of  the  corresponding  linear  problem  (Ax  = b) , or  apply  the  Young's  method 
(see  VARGA  T83],  YOUNG  [84]).  Actually  the  last  method  has  given  very  good 
results  even  for  matrices  A which  are  not  M-matrices. 


Remark  3 . 1 : In  practice,  when  using  S.O.R.  with  truncation  to  solve  dis- 
crete obstacle  problems,  the  first  components  of  xn  to  converge  are  those 
for  which  the  discrete  solution  coincides  with  ip  . It  appears  in  fact  tnat 
most  of  the  computational  time  is  used  to  compute  the  approximate  solution 

of  - Au  = f on  51+  = (x  e 51,  u(x)  > iKx)},  with  u|  ^ = \1> I as  boundary 

'352+  351+ 

conditions.  And  indeed  we  have  observed  that  the  optimal  value  of  w corres- 
ponds to  the  optimal  choice  for  the  approximate  solution  of  the  correspon- 
ding linear  Dirichlet  problem  on  51+. 


Concluding  comments  : The  main  advantages  of  S.O.R.  methods  with  truncation 
arc  that  : 

- they  are  easy  to  program 

- they  require  few  computer  storage, 

- they  have  however  some  drawbacks  which  are  that  they  are  mainly 

limited  to  second  order  potential  problems.  Indeed  they  usually 

show  a fairly  slow  convergence  when  applied  to  the  solution 

2 

of  obstacle  problems  related  to  A or  to  elasticity  operators 
(the  block  variants  of  algorithm  (3 . 3— (3.5)  are  not  so  easy  to 
program  since  the  truncation  has  to  be  replaced  by  a more  compli- 
cated projection  step). 


3.1.3.  Solution  of  the  obstacle  problem  by  duality  methods. 

Several  dual  problems  may  be  associated  to  the  obstacle  problem  (2.1), (2. 2) 
(see  Sec.  4 below  for  one  of  them).  Among  these  dual  formulations  the  follo- 
wing is  well  suited  for  computations. 


Let  define 

2 

A = {y  £ L (51),  y > 0 a.e.}  , 
Vg  = (v  e H*  (51)  , v=g  on  351} 
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and  a lagrangian  functional  sL  : H (fl)  x L (fl)  + E by 


j£(v,u)  = j 


|Vv|2dx- 


fv  dx  - y(v~ip)dx  . 

Q h 


Assume  that  the  solution  u of  (2.1), (2. 2)  belongs  to  H (fl)  , then  X = -Au-f 
is  the  unique  solution  of  the  dual  problem 


(3.6)  Max  Min  j£(v,y)  ; 

ye  A veV 

8 

moreover  we  can  easily  prove  that  {u,A}  is  the  unique  saddle-point  of  tt  over 

V xA  . 
g 


From  these  properties  we  can  use  the  Uzawa's  algorithm  discussed  in  GLOWINSKI- 
LIONS-TREMOLIERES  [ 57,  Ch.  2)  which  takes  the  following  form  : 

(3.7)  A°  e L 2(fl),  arbitrary  given  (X°  = 0 for  ex.) 

then  for  n>0,  assuming  An  known  we  obtain  un  and  An+'  by 

i-  Aun+I  = f + An  in  fi  , 
uI1+1  = g on  3(1, 

(3.9)  An+1  = P+(An+p(t|j-un) ) , p>0  , 


where 


P+(y)  = Sup(0,y) 

Using  GLOWINSKI-LIONS-TREMOLIERES  [57,  Ch.  2]  we  can  prove 

Theorem  3.2  : Assume  that  u e (fi)  , then  V A°  e ((1)  we  have 

lim  |(  u11— u 1 1 = 0 

n -*■  + °°  H (d) 

if  0 < p < 2 g 
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is  the  smallest  eigenvalue  in 


where  £ 


l - Aw  = f w _iri  SI , 
' w = 0 on  3fi . 


Comment  s : 

1)  If  one  uses  A°  = 0 we  have  observed  that  the  smallest  is  = {x  e f),  u(x)  =l(/(x) } , 
the  fastest  is  the  convergence. 


2)  The  above  method  (in  its  discrete  form)  is  very  well  suited  to  users  having 

at  their  disposal  finite  difference  or  finite  element  elliptic  solvers. 




3)  Another  advantage  of  this  duality  approach  is  that  it  gives  directly  A 
whose  Mechanical  interpretation  is  interesting  since  it  is  the  reaction 
force  of  the  obstacle  on  the  membrane. 

4)  Valiants  of  the  above  algorithm  have  been  successfully  used  witii  A replaced 
by  A"  or  elasticity  operators. 


3.1.4.  Solution  of  the  obstacle  problem  by  penalty  methods. 

Like  in  Sec.  3.1.3.  we  shall  focus  our  attention  to  the  continuous  obstacle 
problem  whose  formalism  is  simpler. 


With  V like  in  Sec.  3.1.3.  we  consider 

8 


(3.10)  Min 

veV 


Vv  dx  - 


fv  dx  + j (x) } 


8 


with 


jC(v)  = 2i 


(3.11) 

where  q + = sup  (0,q) 


(tj;-v)  + | 2dx 


The  minimization  problem  (3.10)  is  in  fact  equivalent  to  the  Non  linear  Dirichlet 

p roblem 


i - Auf  - ~ )+  = f in  Q, 

(3.12) 

! n(  = g on  311. 

Concerning  the  convergence  of  u^  to  the  solution  u of  (2.1)  (2.2)  we  can 
pr  ive  that 

lim  ||  u -u  ||  = 0 . 

£+0  b H (fi) 

The  non  linear  problem  (3.10),  (3.12)  (in  fact  its  discrete  variants)  can 
be  solved  by  various  methods,  like  e.g.  Non  Linear  S.O.R. .Conjugate  Gradient 
with  Scaling  (cf.  CONCUS-GOLUb-O’ LEARY  [37],  DOUGLAS-DUPONT  [ 44J)  using  as 
scaling  operator  a discrete  form  of  -A,  or  an  operator  obtained  by  Incomplete 
Cholesky  Recomposition  (see  ME1RINJK-VAN  DE  VORST  [ 68 J for  details). 

Comments  : The  main  inconvenient  of  that  penalty  method  is  that  it  requires 
a small  e to  have  u^-u  small.  Then  (3. 1 0) , (3. 1 2)  are  ill-conditioned  problems. 
However  this  technique  has  been  successfully  used  in  Optimal  Control  or 
Optimal  Design  problems  in  which  the  state  equation  is  replaced  by  a varia- 
tional inequality. 

3.1.3.  Other  methods. 

The  above  obstacle  problem  can  also  be  solved  by  Conjugate  gradient  with 
truncation  methods  (see  e.g.  TREMOLIERES  [81])  or  by  using  a discrete  time 
dependent  approach  requiring  at  each  time  step  the  solution  of  a problem 
of  similar  type  but  better  conditioned. 

We  can  also  use  augmented  lagrangian  methods  (in  this  direction  we  refer  to 
CHAN-GLOWINSKI  [32]). 

3.2.  Iterative  solution  of  the  elasto-plastic  torsion  problem. 

In  this  section  we  just  describe  one  algorithm  and  send  the  reader  to  GL0WINS1 
LIONS-TREMOLIERES  [57,  Ch.  3],  CEA-GLOWINSKI-NEDELEC  [3|],  FORTIN-GLOV INSKI 
[51  , Ch.  3]  and  GL0WINSKI-MARR0CC0  [58]  for  more  details  and/or  other  methods 

The  problem  under  consideration  is  (2.29)  ; we  associate  to  this  problem  a 

1 °o  _ , 

lagrangian  functional  £ : H^(ft)  *L  (ft)  ->  K defined  by 
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* 


(3.13) 


a?(v,M) 


_l_ 

9 


Vv| ~dx 


fv  dx  + | u(|Vv|"-l)dx  . 

'fl 


OO  -1  . . , 

If  A = {y  e L (Q)  , \i  ^ 0 a.e.  then  we  can  easily  prove  that  n lu,A,  is  a 
saddle-point  of  ^ over  H^(il)x A,  then  a is  the  solution  of  (2.29).  If  f = const.. 

the  existence  of  such  a saddle-point  has  been  proved  by  H.  BREZTS  i 16 j,  for 
more  general  f the  situation  is  not  clear  at  the  moment.  However  the  existence 
of  such  saddle-point  is  proved  for  the  finite  difference  or  finite  element 
approximation  of  the  elasto-plastic  problem  (see  GL0WINSK1- LIONS-TREMOLIERES  157, 
Ch.  3]  and  CEA-GLOWINSKI-NEDELEC  [31]  for  further  information).  From  the 
above  properties  it  is  then  natural  to  use  the  following  algorithm  to 
compute  the  saddle-points  of  s£  and  therefore  u : 


(3.14)  X°  £ A , arbitrarily  given  , 


then  for 


n s 0 compute  u and  X 


n+ 1 


by 


(3.15) 


(3.16) 


- V* (l+Xn)Vun  = f in  0 , 


Xn+'  = P+(Xn  + p(|Vun|2-l))  , p >0  , 


with  P+  like  in  Sec.  3.1.3. 

For  the  convergence  of  (3 . 1 4) -(3 . 1 6)  see  the  two  references  above. 

We  observe  that  the  solution  of  (3.15)  require  at  each  iteration  the  solution 
of  a Dirichlet  problem  whose  right-hand  side  depends  upon  n (via  Xn) . However 
(3. 14)— (3. 16)  appears  as  an  efficient  algorithm  for  solving  (2.29). 

3.3.  Iterative  solution  of  the  Bingham  flow  problem  in  a pipe. 

We  consider  in  this  section  the  iterative  solution  of  (2 . 38) , ( 2 . 39)  . Actually 
the  simplest  method  to  solve  this  problem  is  based  on  the  characterization 
(2.41 ),  (2.42)  . The  algorithm  is  the  following 

(3.17)  p°  € L2(0)  xL2(fi)  is  arbitrarily  given  (,pf  = 0 for  ex.) 


then  for  n-  0 one  defines  un  and  Xn  ' from  Xn  by 
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(3. 18) 


. n _ _ n . 

- ViAu  = f + gv*p  in  £2, 

= 0 ’ 

(3.19)  pn+1  = PA  (pn+  pVu11)  , p > 0 , 

where  A = { q e (ft)  x (ft)  , | q (x)  | < 1 a. e . } and 

V>  ' ,up(T, up  v,«i.2(fi)xi.2<n)  . 

It  follows  then  from  CEA-GLOWINSKI  [30]  and  GLOWINSKI-LIONS-TREMOLIERES  [57, 
Ch.  2 and  5] 


Theorem  3.3  : Assume  that  we  have 

0 < p<  - 

g 

then  V p°  £ L (ft)  x L (ft)  we  have 

lim  ||un-u||  = 0 , 

H (ft) 

n oo  oo 

lim  p = p _in  L (ft)  x l (ft)  weak-*  , 
n -*■  +°° 

where  u is  the  solution  of  (2.38) , (2.39)  and  where  p is  such  that  {u,p}  obeys 
(2.41),(2.42). 


Remark  3.2  : The  above  function  p is'  actually  a solution  of  the  dual  problem 


Max  Min  { y 
qeA  velr(ft) 


ft 


Vv|  dx  - | fv  dx  + g 

ft 


ft 


q*Vv  dx)  . 


Remark  3.3  : More  efficient  algorithms  based  on  augmented  lagrangian  techniques 
are  described  in  GABAY-MERCIER  [52],  FORTIN  [50],  FORTIN-GLOWINSKI  [51], 
GLOWINSKI-MARROCCO  [58].  These  algorithms  are  more  complicated  to  handle  than 
(3 . 1 7) -(3 . 1 9)  which  is  definitely  the  simpler  efficient  method  to  solve  (2.38), 
(2.39)  and  their  discrete  variants. 
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•t.  - A 1'PROX  1MAT i ON  OF  THE  OBSTACLE  PROBLEM  BY  MIXED  FINITE  ELEMENT  MKTHODS . 

A . I . Orientation. 

Following  BREZZI-HAGER-RAV1ART  [22]  we  shall  consider  in  this  section  mixed 
finite  element  approximations  of  the  obstacle  problem  (2.1),  (2. 2)  (see  also 
for  related  studies  HI.AVACHEK  161]). 

For  more  complicated  variational  inequalities,  solved  also  by  mixed  finite 
element  methods  we  refer  to,  e.g.,  JOHNSON  [62],  [63],  BREZZI-J0HNS0N-MERC1ER 

[23]  , BEGIS-GLOWINSKI  [8]. 

A. 2.  A dual  formulation  of  the  obstacle  problem. 

It  is  clear  that  (2.1), (2. 2)  has  also  the  following  formulation 

(A.l)  Min  { y |q|"dx  - [ fv  dx} 

1 v , q } e K 'ft 

where 

2 N 

(4-2)  K = { { v , q}  € H (ft) x ( L (ft) ) W , q=Vv  , v=g  on  3ft,  v S a.e.  on  ft). 


We  associate  to  (4.1), (A. 2)  the  lagrangian 

a£(v,q,y)  = y 


( |<l|2dx  - 

fv  dx  + 

h J 

ft 

y • (Vv-q) dx 
ft 


and  then  consider  the  dual  problem  of  (2.1), (2. 2)  related  toad  , i.e. 

(A. 3) 


where 


Max  Min  ad(v,q,y) 

ye(I.2(ft))N  {v.qkft 


K = • ■-  v,q}  e H (ft)  x L^(ft)  , v=g  on  3ft,  v > t|/  a.e . on  ft} 


Fortunately  the  explicit  form  of  (A. 3)  is  known  and  is  as  follows 


(A. A) 


with 


-y  1 i q ! ^dx  + V« 

qeC  Jft  'ft 


q dx  - 


g q *n  df) 
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(4.5) 


C = {q  t H(div,12)  , V*q  + fs  0 a.e.  on  H) 


where 


H(div.n)  = { q £ (L2(f2))N  , V-qcL2®} 


Conversely  (4.4)  has  a unique  solution  p such  that  p = Vu  where  u is 
solution  of  the  obstacle  problem  (2.1),  (2.2). 

Def ine 


A = {y  e L (11)  , y>0  a.e.  on  11} 


then 


(4.6) 


q e c q £ H(div ,12)  , 


(V*q+f)y  dx  < 0 V y £ A 


n 


4.3.  A mixed  approximation  of  the  dual  problem  (4. 4), (4. 5) 

2 

We  still  assume  that  H is  a polygonal  domain  of  R and  SC=>.  a standard 

h 

lation  of  H . Let  consider  the  following  approximations  of  H(div,li), 
A,  C respectively 

2 2 

= {qh  e L (H)  x L (H)  , qh*n| 9T  e sidew  ise  where  n is 

unit  normal  along  3T,  V*qh|TePk} 

(from  the  above  properties  H,  c H(div  ,11) ) , 

L„  - (|lhiL!«)  , UhlT  € Pk  VTtffh)  , 

A^  = A n L^  = { Uv,  e , ^>0  a.e.  on  H}  , 

(V*qVi+f)y>i  dx  < 0 Vy^eA^} 


Ch  = {qheHh  * 


An  approximate  problem  of  (4.4),  (4.5)  is 


the 


triangu- 
L2  (11)  , 


Liu 


-42- 


-43- 


5.  - FURTHER  COMMENTS.  CONCLUSION 

The  variational  inequality  methodology  appears  as  an  efficient  tool  for 
solving  problems  which  may  be  originally  formulated  in  a more  classical 
way.  Let  us  consider  two  examples  of  such  a situation. 

Example  1 - A family  of  mildly  non  linear  elliptic  problems. 

N . 1-1 

Let  SI  be  a bounded  domain  of  R with  smooth  boundary,  A : H (ft)  H (ft) 

1 | Q 

(=  (HQ(ft))')  a strongly  elliptic  isomorphism,  feH  (ft)  , $ : R -►  R , 4>eC°, 

non  decreasing  (we  can  always  suppose  that  <|>  (o)  =0).  We  consider  then  the 

following  non  linear  elliptic  problem  (of  monotone  type) 


(5.1) 


Au  + <jj(u)  = f . 


Let  define  a : H*  (ft)  xH*(ft)  -*  R , bilinear  continuous  and  H*  (ft) -elliptic  by 
o o o 


a(v,w)  = <Av,w> 


where  <*,•>  denote  the  duality  pairing  between  H '(ft)  and  h' (ft) . Let  define 

1 ° 
also  j : H (ft)  + R by 


j (v)  = <I>(v)dx 


where  $(t)  = <f>(x)dT  • It  is  clear  that  j(*)  is  convex,  proper,  l.s.c. 

I o 

Actually  in  view  of  solving  (5.1)  it  is  convenient  to  consider  the  following 
(E.V.I.). 


(5.2) 


Find  ueH  (ft)  such  that 
o 


a(u,v-u)+j  (v)-j  (u)  > <f,v-u>  VveHo(ft). 


In  CHAN-GLOWINSKI  [32]  (see  also  GLOWINSKI  [55])  one  proves  that  (5.2)  has 

a unique  solution  which  is  also  the  unique  solution  of  (5.1)  in  H1 (ft) . We 
o 

consider  also  in  the  two  above  references  the  approximation  of  (5.1), (5. 2) 
by  piecewise  linear  finite  elements.  This  is  precisely  a situation  in 
which  the  convergence  results  of  Sec.  1.6.1.  cannot  be  applied  (directly  at 
least),  since  j(*)  is  not  continuous  on  H^(ft)  in  general.  However  the  strong 
convergence  of  the  approximate  solution  can  be  proved  (see  the  above  references 
for  more  details). 


J 


Example  2 - Transonic,  potential  flows. 

This  problem  is  far  more  important  and  difficult  than  the  problem  of  Example  1. 
In  fact  it  is  not  an  elliptic  problem  unless  the  flow  remains  purely  subsonic. 
The  problem  is  to  find  a velocity  potential  <p  defined  on  (2  (the  domain  of 
the  flow)  such  that 


(5.3) 


| V*  (p (<t>) V4>)  = 0 in  SI, 

( + convenient  boundary  conditions, 


with 


pW)  o . JaLp/vi 


14  c2 

y- 1 * 


where  p = const.,  y = 1 .4  in  air  , = critical  velocity.  The  flow  velocity 

is  given  by  v = Vf}'. 


If  fl  is  not  simply  connected,  <|>  has  to  obey  a circulation  condition  given  by 
the  Kutta-Joukowsky  condition  (see  LANDAU-LIFCHITZ  [64]  for  more  details,  and 
also  the  ref erences  be  low) . Actually  the  above  relations  are  not  sufficient  to 
obtain  only  physical  solutions,  i.e.  solutions  without  expansion  shocks  ; to 
avoid  these  non  physical  solutions  an  Entropy  Condition  has  to  be  prescribed. 
We  have  found  convenient  to  require 

(3.4)  (A<}>)+  e LP(S1)  Vp>l  . 


The  more  common  values  of  p are  p=2  and  p=  + 00 . In  order  to  solve  (5.3),  taking 

(5.4)  into  account  we  have  introduced  the  variational  problem 


(5.5) 


Min  { 
w«  X 


| V4>(w) 
ft 


"dx  + y 


(Aw-C) + | ^dx} 


where  X is  a set  (usually  convex)  of  admissible  velocity  potentials,  y a posi- 
tive parameter,  C a given  constant  (or  function)  and  <J>(w)  the  solution  of  the 
elliptic  problem 


(3.6) 


A<$(w)  = V*  (p(w)Vw)  in  ft, 


+ boundary  conditions  + (possibly)  Kutta-Joukowsky  conditions. 

In  fact  (5. 5), (5. 6)  is  a least  square  formulation  of  (3.3)  taking  ( 3 . 4 > 
into  account.  The  problem  (5.5)  which  is  indeed  a non  linear  fourth  ordc 
variational  problem  has  been  solved  by  a mixed  finite  element  method  couple 
to  a conjugate  gradient  algorithm  with  scaling.  For  further  details  we  refej 
to  GLOWINSKI-PIRONNEAU  [59],  [60],  BRISTEAU  [24],  [ 25  J , BRISTEAU-GL0W1NSK1 
PERRIER-PERIAUX-PIRONNEAU-POIRIER  [26]  where  references  to  other  methods 
for  solving  (5.3)  are  also  given. 

To  conclude  this  survey  on  Elliptic  Variational  Inequalities,  let  us  mt  o 
several  books  or  reports  relevant  to  the  subject  : 

— DUVAUT- LIONS  [45]  fur  the:  mathematical  aspect  and  applicat  o * 
to  Mechanics  and  Physics. 

- CEA  [28],  GLOWINSKI-LIONS-TREMOLIERES  [57],  GL0WINSK1  >531, 
[54]  , [55],  TREMOLIERES  [81 J for  the  numerical  analysis 
(approximation  and  iterative  solution) . 

- BENSOUSSAN-LIONS  [9  ] for  Quasi  Variational  Inequalities 

- BAIOCCHI-CAPELO  [6  ],  for  ^pplicati  ons  to  the  solution  of 
the  free  boundary  problems  related  to  flows  in  porous  meu  j . 
(for  that  last  subject  see  also  BAIOCCHI-COMINCIOLI-MAGLM  : 
POZZI  [7],  CRYER-FETTER,  loc.  cit.). 
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